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Abstract. We present an extraction of the E2/M1 ratio of the Z\(1232) from experimental data applying 
an effective Lagrangian model. We compare the result obtained with different nucleonic models and we 
reconcile the experimental results with the Lattice QCD calculations. 

PACS. 14.20.Gk Baryon resonances with S = - 25.20.Lj Photoproduction reactions - 13.60.Le Meson 
production 



The deformation of the nucleon and its first excita- 
tion, the .4(1232), is a topic that has focused the atten- 
tion of many researchers in the last years from both the 
experimental and the theoretical sides [T]. The possibility 
of such deformation has been studied using the E2/M1 
ratio (EMR) of the jN _> 4(1232) transition [2J. The 
emission (absorption) of a photon by a spin-3/2 particle 
' involves a magnetic dipole (Ml) multipolarity and an elec- 
tric quadrupole (E2) multipolarity. From experiments it is 
found that E2 is small but not zero, which evokes a de- 
formed nucleon picture. A deviation from zero of the EMR 
is a clear indication of the existence of such deformation 
and allows to quantify it. This ratio is mainly obtained in 
two different ways, from nucleonic models such as quark 
models or Lattice QCD, and from experimental data. Rec- 
onciliation of both extractions is significant in order to un- 
derstand the structure of the nucleon. This work is con- 
cerned with the extraction of the intrinsic E2/M1 ratio 
of the 4(1232) from experimental data using a reaction 
model. 

In [3l|4] we have developed a pion photoproduction 
model up to 1 GeV of photon energy based upon effective 
Lagrangians. The model follows closely the work of Gar- 
cilazo and Moya de Guerra [5] and has similarities with 
the work of Sato and Lee [B] as well as with other works 
[7HBUH] based on the seminal work of Peccei [TU] ■ The re- 
action model allows us to isolate the contribution of the 
4(1232) - taking into account the high energy behaviour 
of the tail of the resonance -, to calculate its EMR, and 
to compare it with the values provided by nucleonic mod- 
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els. In addition to Born terms (those which involve only 
photons, nucleons, and pions) and vector meson exchange 
terms (p and lu exchanges), the model includes all the 
four star resonances in Particle Data Group (PDG) [19] 
up to 1.7 GeV mass and up to spin-3/2: 4(1232), N(1440) J , 
N(1520), 4(1620), N(1650), and 4(1700). 

The model displays chiral symmetry, gauge invariance, 
and crossing symmetry, as well as a consistent treatment 
of the interaction with spin-3/2 particles that avoids well- 
known pathologies present in previous models OH]. The 
dressing of the resonances is considered by means of a phe- 
nomenological width which takes into account decays into 
one 7r, one 77, and two tt. The width fulfills crossing sym- 
metry and contributes to both direct and crossed channels 
of the resonances. 

We assume that the final state interactions (FSI) fac- 
torize (nN rescattering) and can be included through the 
distortion of the irN final state wave function. The calcula- 
tion of the distortion requires one to calculate higher order 
pion loops or to develop a phenomenological potential FSI 
model. Both approaches are far of the one we apply. The 
first one is overwhelmingly complex and the second would 
introduce additional model-dependencies, which are to be 
avoided in the present analysis, because we are mainly 
interested in the bare properties of the resonances. We 
rather include FSI in a phenomenological way by adding 
a phase <5fsi to the electromagnetic multipolcs. We de- 
termine this phase so that the total phase of the electro- 
magnetic multipole is identical to the one of the energy 
dependent solution of SAID [11] . In this way we are able 
to isolate the electromagnetic vertex and remove the FSI 
effects. 
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Fig. 1. M 1+ and E 1+ electromagnetic multipoles. Curve conventions: solid, real part of the multipole; dashed, imaginary part 
of the multipole. Data taken from [1 1 J . 



In order to assess the parameters of the model we had 
to minimize the function \ 2 defined by 



E 



Mf (Ai 



AMf p 



where M. exv stands for the current energy independent 
extraction of the multipole analysis of SAID up to 1 GeV 



7T°p 



process (TT]. AA4 exp is the error and Ai th is the multi- 
pole given by the model which depends on the parameters 
Ai, . . . , A„, which stand for the electromagnetic coupling 
constans of the resonances and the cutoff A which regu- 
larizes the high energy behaviour of the Born terms. The 
masses and the widths of the resonances have been taken 
from the multichannel analysis of Vrana, Dytman, and 
Lee [H] and from a speed plot calculation [3] . The EMR 
we present has been obtained as an average of the results 
obtained with both sets of masses and widths. Electro- 
magnetic multipoles are complex quantities and we have 
taken into account 763 data for the real part of the multi- 
poles and the same amount for the imaginary part. Thus, 
m = 1526 data points have been used in the fits. 

In order to fit the data and determine the best pa- 
rameters of the resonances we have written a genetic al- 
gorithm combined with the E04FCF routine from NAG 
libraries [13]. Although genetic algorithms are computa- 
tionally more expensive than other algorithms, in a min- 
imisation problem it is much less likely for them to get 
stuck at local minima than for other methods, namely 
gradient based minimisation methods. Thus, in a mul- 
tiparameter minimisation like the one we face here it is 
probably the best possibility to search for the minimum 
[HUH]. In Fig. [T] we show the fits to Mf{ 2 
tromagnetic multipoles. 



and elec- 



Different definitions of the EMR have been employed 
in the literature. We should distinguish between the in- 
trinsic or bare EMR of the Zl(1232) and the directly mea- 
sured value in experiments which is often called physi- 
cal or dressed EMR value [TI pl IB]. The physical EMR 
is obtained as the ratio between the imaginary parts of 
eV2 and Aff^ 2 electromagnetic multipoles at the invari- 



for E 0+ , Mi_, E 1+ , M 1+ , E 2 -, and M 2 _ multipoles in 
the three isospin channels / = n for the 7p 
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Since all the reaction models are fitted to the 

experimental electromagnetic multipoles, they generally 
reproduce the physical EMR value within the error bars 
(see tabled]). We obtain 



EMR physical = 



Im 


" p 3/2" 


Im 





x 100% = (-3.9 ±1.1)%. (2) 



However, this measured EMR value is not directly avail- 
able from theoretical models of the nucleon and its reso- 
nances. Instead, if we want to compare to models of nu- 
cleonic structure, it is necessary to extract the bare EMR 
value of Z\(1232) which is defined as 



EMR 



G 



Zi(1232) 
E 



G 



Zl(1232) 
M 



x 100% = (-1.30 ±0.52)%, (3) 



The EMR defined in this way depends only on the 
intrinsic characteristics of the Z\(1232) and can thus be 
compared directly to predictions from nucleonic models. It 
is not, however, directly measurable but must be inferred 
(in a model dependent way) from reaction models. 

The intrinsic quadrupole deformation of the Z\(1232) 
is found to be EMR= (-1.30 ± 0.52) %, indicative of a 
small oblate deformation. In Tables [1] and [2] we compare 
our EMR values (bare and physical) to the ones extracted 
by other authors using other models for pion photopro- 
duction, as well as to predictions of nucleonic models. 
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Table 1. Comparison of EMR physical values from experiments 
compared to the values obtained with reaction models. 

EMR physical Ref. 



Experiments 








LEGS Collaboration 


(" 


-3.07 ± 0.26 ± 0.24) % 


US] 


Al Collaboration 


(" 


-2.28 ±0.29 ±0.20)% 


17 


A2 Collaboration 


(" 


-2.74 ± 0.03 ± 0.30) % 


18 


Particle Data Group 


(" 


-2.5 ±0.5)% 


PS 


Reaction models 








Fernandez-Ramirez et al. 


(' 


-3.9 ±1.1)% 


na 


Pascalutsa and Tjon 


(' 


-2.4 ±0.1)% 


m 


Sato and Lee 




-2.7% 


m 


Fuda and Alharbi 




2.09% 


20 



Table 2. Comparison of EMR barc values extracted from ex- 
periments through reaction models compared to the values ob- 
tained with nucleonic models. 





EMR baro 


Ref. 


Reaction models 






Fernandez-Ramirez et al. 


(-1.30 ±0.52)% 




Pascalutsa and Tjon 


(3.8 ± 1.6)% 


m 


Sato and Lee 


-1.3% 


A 


Davidson et al. 


-1.45% 


7 


Garcilazo and Moya de Guerra 


-1.42% 


A 


Vanderhaeghen et al. 


-1.43% 


9 


Nucleonic models 






Non-relativistic quark model 


0% 


HO 


Constituent quark model 


-3.5% 


22 


Skyrme model 


(-3.5 ± 1.5)% 


M 


Lattice QCD (Leinweber et al.) 


(3 ±8)% 


M 


Lattice QCD (Alexandrou et al.) 

(Q 2 = 0.1 GeV 2 , m n = 0) 

(Q 2 = 0.1 GeV 2 , = 370 MeV) 




m 


(-1.93 ±0.94)% 




(-1.40 ±0.60)% 





Reconciliation of the experimental value of the E2/M1 
jN -> zi(1232) transition ratio (EMRP h y sical ) with the 
one obtained using Lattice QCD (EMR barc ) within a con- 
sistent and sound framework (besides our analysis only in 
[8] a consistent treatment of the spin-3/2 interaction is 
performed), is one of the goals of this work [15]. Our re- 
sults also indicate that quark models need improvements 
in order to reproduce the value obtained from experiment. 
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FIS2005-00640 and BFM2003-04147-C02-01. 



References 

1. B. Krusche and S. Schadmand, Prog. Part. Nucl. Phys. 51 
(2003) 399. 

2. A.M. Bernstein, Eur. Phys. J. A 17 (2003) 349. 

3. C. Fernandez- Ramirez, Electromagnetic production of light 
mesons. PhD dissertation. Universidad Complutense de 
Madrid (2006). 



4. C. Fernandez- Ramirez, E. Moya de Guerra, and J.M. Udi'as, 
Ann. Phys. (N.Y.) 321 (2006) 1408. 

5. H. Garcilazo and E. Moya de Guerra, Nucl. Phys. A 562 
(1993) 521. 

6. T. Sato and T.-S.H. Lee, Phys. Rev. C 63 (2001) 055201 

7. R.M. Davidson, N.C. Mukhopadhyay, and R.S. Wittman, 
(1991) Phys. Rev. D 43 (1991) 71. 

8. V. Pascalutsa and J.A. Tjon, Phys. Rev. C 70 (2004) 
035209. 

9. M. Vanderhaeghen, K. Heyde, J. Ryckebusch, and M. Waro- 
quier, Nucl. Phys. A 595 (1995) 219. 

10. R.D. Peccei, Phys. Rev. 181, 1902 (1969). 

11. R.A. Arndt, W.J. Briscoe, I.I. Strakovsky, and R.L. Work- 
man, Phys. Rev. C 66 (2002) 055213, SAID database, 
http : //gwdac . phys . gwu . edu 

12. T.P. Vrana, S.A. Dytman, and T.-S.H. Lee, Phys. Rep. 
328 (2000) 181. 

13. Numerical Algorithms Group Ltd., Wilkinson 
House, Jordan Hill Road, Oxford OX2-8DR, UK. 
http : //www . nag . co . uk/ 

14. D.G. Ireland, S. Janssen, and J. Ryckebusch, Nucl. Phys. 
A 740 (2004) 147. 

15. C. Fernandez-Ramirez, E. Moya de Guerra, and J.M. 
Udi'as, Phys. Rev. C 73 (2006) 042201 (R). 

16. G. Blanpied et al, Phys. Rev. C 64 (2001) 025203. 

17. S. Stave et al., arXiv:micl-ex/0604013 

18. J. Ahrens et al., Eur. Phys. J. A 21 (2004) 323. 

19. W.-M. Yao et al., J. Phys. G 33 (2006) 1, 
http: //pdg. lbl .gov 

20. M.G. Fuda and H. Alharbi, Phys. Rev. C 68 (2003) 064002. 

21. CM. Becchi and G. Morpurgo, Phys. Lett. 17 (1965) 352. 

22. A.J. Buchmann, E. Hernandez, and A. Faessler, Phys. Rev. 
C 55 (1997) 448. 

23. A. Wirzba and W. Weise, Phys. Lett. B 188 (1987) 6. 

24. D.B. Leinweber, T. Draper, and R.M. Woloshyn, Phys. 
Rev. D 48 (1993) 2230. 

25. C. Alexandrou, Ph. de Forcrand, H. Neff, J.W. Negele, 
W. Schroers, and A. Tsapalis, Phys. Rev. Lett. 94 (2005) 
021601. 



